Compact stars within an asy-soft quark-meson-coupling model 
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We investigate compact star properties within the quark meson coupling model (QMC) with a 
soft symmetry energy density dependence at large densities. In particular, the hyperon content 
and the mass/radius curves for the families of stars obtained within the model are discussed. The 
hyperon-meson couplings are chosen according to experimental values of the hyperon nuclear matter 
potentials, and possible uncertainties are considered. It is shown that a softer symmetry energy gives 
rise to stars with less hyperons, smaller radii and larger masses. Hyperon-meson couplings may also 
have a strong effect on the mass of the star. 
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I. INTRODUCTION 

In the last years important efforts have been done to 
determine the density dependence of the symmetry en- 
ergy of asymmetric nuclear matter (see the reviews [IHl] 
and references therein). Correlations between different 
quantities in bulk matter and finite nuclei have been es- 
tablished. For instance, the correlation between the slope 
of the pressure of neutron matter at p = 0.1 fm -3 and 
the neutron skin thickness of 208 Pb [3, [||, or the corre- 
lation between the crust-core transition density and the 
neutron skin thickness of 208 Pb @ are well determined. 
Presently, there also exist different experimental mea- 
surements that constrain the saturation properties of the 
symmetry energy [7J. 

The quark-meson-coupling (QMC) model [8l-flfj| is an 
effective nuclear model that takes into account the inter- 
nal structure of the nucleon explicitly. Within the QMC 
model, matter at low densities and temperatures is a sys- 
tem of nucleons interacting through meson fields, with 
quarks and gluons confined within MIT bags [11]. For- 
matter at very high density or temperature, one expects 
that baryons and mesons dissolve and that the entire sys- 
tem of quarks and gluons becomes confined within a sin- 
gle, big, MIT bag. Within QMC it is possible to describe 
in a consistent way both nucleons and hyperons [l2| . The 
energy of the baryonic MIT bag is identified with the 
mass of the baryon and is obtained self-consistently from 
the calculation. It is important to stress that within the 
QMC model the coupling of hyperons to the cr-meson is 
fixed at the level of the saturation properties of the equa- 
tion of state (EOS). Hypernuclei properties [T3T - [l6| will 
then allow the determination of the coupling of hyperons 
to the isoscalar- vector meson without any ambiguity ex- 
cept for the uncertainty on the experimental hypernuclei 
data. Within the non-linear Walecka models (NLWM) 
this is not possible and some other constraint must be 
imposed, such as using the SU(6) symmetry to fix the 



hyperon- vector meson couplings [17] . 

Recently, new data on neutron stars have been ob- 
tained 0, [Til that theory should explain, namely, the 
large value 1.97 ± 0.04 M of the recent mass measure- 
ment of the binary millisecond pulsar PSR J16142230 
18], and the empirical EOS obtained by Steiner et al. 
from a heterogeneous set of seven neutron stars with 
well-determined distances [l|| which predicts quite small 
radii. These last results, however, should still be con- 
sidered with care because there are many uncertainties 
involved. 

The symmetry energy at saturation is quite well estab- 
lished, however, the density dependence of the symmetry 
energy is not so well known and different models predict 
a wide range of values for the symmetry energy slope at 
saturation. Although the symmetry energy slope of QMC 
at saturation (94 MeV [2(|) is within the range of values 
compatible with experimental observations [211 ] . most of 
the experimental observables that have been proposed 
to obtain a measure of the symmetry energy slope, pre- 
dict smaller slopes Q , which could be as low as 30 MeV. 
Moreover, the large value of the QMC slope prohibits 
the prediction of small radii as the ones indicated by the 
empirical EOS [H. 

In [2(| we have introduced the <5-meson in the QMC 
model and have studied its effect on the density depen- 
dence of the symmetry energy and subsaturation instabil- 
ities of nuclear matter. However, the introduction of the 
(5-meson gives rise to a stiffcr symmetry energy and this 
mechanism will not allow us to obtain a softer symmetry 
energy for the QMC. 

In the present study we consider an extension of the 
QMC that includes a nonlinear term involving the uj and 
p mesons. This term affects the isovector channel of the 
QMC equation of state, namely the density dependence 
of the symmetry energy [6|,[22|, and choosing the coupling 
constant adequately it is possible to correct the stiff be- 
havior of the symmetry energy at large densities in the 
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QMC model. 

Stellar matter within the modified QMC model will be 
studied. We expect that smaller values of the symmetry 
energy slope will give rise to smaller star radii: this has 
been shown both for nucleonic stars [f| [H[ and for hy- 
peronic stars [24]]. In particular, we want to investigate 
how the hyperon content of a compact star is affected 
by the density dependence of the symmetry energy. This 
was recently done within the NLWM including the non- 
linear up term [24J and it was shown that the hyperon 
content could be affected by the density dependence of 
the symmetry energy, and that the radius of stars with 
a mass 1-1.4 Mq increases linearly with the slope, while 
the mass is not affected. However, the QMC model gives 
rise to a softer EOS and generally predicts smaller hy- 
peron fractions in stellar matter. This could affect the 
behavior of the star properties and its relation with the 
slope of the symmetry energy. 

The paper is organized as follows: in section II an ex- 
tension of the QMC model to include the ui — p coupling 
is discussed, in section III results are presented and dis- 
cussed and the final conclusions are drawn in the last 
section. 



II. THE QUARK-MESON COUPLING MODEL 

In what follows we present a review of the QMC model 
and its generalization to include the isoscalar-isovector 
u) — p coupling. 

In the QMC model, the nucleon in nuclear medium is 
assumed to be a static spherical MIT bag in which quarks 
interact with the scalar (a) and vector (u>, p) fields, and 
those are treated as classical fields in the mean field ap- 
proximation (MFA) [1, Q. The quark field, tf) qi , inside 
the bag then satisfies the equation of motion: 



i p - to* - figl u + g% q b 03 ) j> 9t = 



(1) 



where q — u,d,s, to* — m a q — g q a with m q the current 
quark mass, and g q , and g q p denote the quark- meson 
coupling constants. The energy of a static bag describing 
baryon i consisting of three quarks in ground state is 
expressed as 
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where il qi = 



2 . + (Ri to*) 2 , Ri is the bag radius of 

baryon i, x qi is the dimensionless quark momentum, Zi 
is a parameter which accounts for zero-point motion of 
baryon i and is the bag constant. The effective mass 
of a baryon bag at rest is taken to be M* = -E, as . The 
equilibrium condition for the bag is obtained by mini- 
mizing the effective mass, M* with respect to the bag 
radius 
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211.30306 MeV and R, = 0.6 



fm. The unknowns Zi are given in [K 



A. QMC with coupled uj — p fields 

The consideration of a coupling between the isoscalar 
and isovector fields is carried out much like in the manner 
it was performed in [||. A note is however on demand: 
we start out from quarks, which find themselves confined 
in a bag [ll|, and the boundary conditions for achieving 
confinement hold the same. The relevant changes (of 
couplings) and fittings (to the symmetry energy) are done 
otherwise for hadronic matter. The total energy density 
of the nuclear matter then reads 
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where t%B is the isospin projection of baryon B. For the 
nucleons we take x^b = XpB = 1- The corresponding co- 
efficients for the hyperons will be discussed later. In the 
above expression for the energy density, we have intro- 
duced the lu — p couplings. The chemical potentials, nec- 
essary to define the /3-equilibrium conditions, arc given 
by 
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In the above expressions the mean fields for mesons are 
determined by the equations 
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^■Kg p gi^Q, and g^ = 3g% and g p = g q p . 

The model parameters arc obtained by fitting the nu- 
cleon mass and enforcing the stability condition for the 
bag in free space. The desired values of En = e/ p — M = 
— 15.7 MeV at saturation p = po = 0.15 fm -3 , are 
achieved by setting g q = 5.981, g u = 8.954. We take the 
standard values for the meson masses, namely m a — 550 
MeV, m u = 783 MeV, and m p = 770 MeV. The cou- 
plings g p and A„ are determined so that £ sym = 23.27 
MeV at p = 0.1 fm~ 3 (k F - 1.14 fm -1 ). The parameters 
g p and A„ are listed in Table HI 



III. RESULTS AND DISCUSSIONS 

Before applying the modified QMC model to the study 
of stellar matter we discuss its properties at saturation 
and subsaturation densities. 
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TABLE I. Nuclear matter properties of the models used in the 
present work. All quantities are taken at saturation, where 
B/A=15.7 MeV and the compressibility K = 290 MeV. p t 
and Y p t are the density and proton fraction at the crust-core 
transition estimated from the thermodynamical spinodal sec- 
tion. 



Model 


A„ 


9p 




L 


Pt 










(MeV) (MeV) (fm" 3 ) 




QMC 





8.8606 


33.70 


93.59 


0.076 


0.022 


QMCojp 


0.01 


8.9837 


33.02 


84.99 


0.079 


0.024 




0.02 


9.1122 


32.43 


77.35 


0.081 


0.025 




0.03 


9.2463 


31.88 


70.55 


0.083 


0.026 




0.05 


9.5335 


30.87 


59.03 


0.089 


0.029 




0.10 


10.3869 


27.78 


39.04 


0.098 


0.033 


NL3 [25] 





8.9480 


37.34 


118.30 


0.065 


0.021 


TW[26] 
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FIG. 1. (Color online) Symmetry energy (a), its slope pa- 
rameter L = 3po£' S ym (b), and curvature K sym for QMC and 
QMCojp for different values of the coupling A„ . 



The saturation properties of nuclear matter for all the 
QMC models considered in the present work are shown 
in Table HI For comparison, we also include the proper- 
ties of two well known relativistic nuclear models NL3 
[25| and TW [Ifij]. Except for the the symmetry energy 
of the QMCwp model with A„ = 0.1, which could be 
a bit too low, both the symmetry energy and its slope 



L = 3pod£ sym / 'dp of all the parametrization are well 
within the experimental constraints coming from differ- 
ent sources p}. 

FigJT] shows the dependence on the density of the 
symmetry energy £ sym and its slope for the different 
parametrizations of the non-linear uj — p term. The sym- 
metry energy is given by 
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Mq 2 and Mq is the nucleon effec- 



where ep — y n> F 
tive mass in symmetric nuclear matter. The symmetry 
energy within QMC shows a rather linear behavior with 
density. This is a feature of many NLWM models. The 
nonlinear up term changes the density dependence of the 
symmetry energy, and, as a result, the symmetry energy 
of the QMCuip (FigOJa)) becomes softer at higher densi- 
ties. This is confirmed by the slope parameter L, plotted 
in Fig. [TJb): above ~ 0.06 fm~ 3 , L becomes smaller the 
larger the coupling A„ is. However, below p = 0.1 fm~ 3 
the symmetry energy is larger in the models with a larger 
A v and this has important effects on the properties of the 
crust-core transition. 



A. Crust-core transition 

The QMCw/9 presents larger instability regions than 
QMC at subsaturation densities and large isospin asym- 
metries. The larger the magnitude of the coupling A„ the 
smaller the slope L and the larger the instability region. 
The same behavior was obtained in_J22| for NL3wp and 
is contrary to the one obtained in [20| with the inclusion 
of the <5-meson in the QMC model. In (2?J the effect of 
the slope L on the spinodal surface at large asymmetries 
was discussed and it was shown that larger values of L 
give rise to smaller spinodal regions at large asymme- 
tries, where matter is closer to neutron matter. Neutron 
matter pressure is essentially proportional to the slope L 
and, therefore, a larger L corresponds to a harder EOS. 
The crust-core transition densities are shown in Table I 
and it is seen that they increase with the increase of the 
lj — p coupling. The proton fraction at the transition den- 
sity also tends to increase QMC gives a higher transition 
density than the one obtained within NL3. For A„ in the 
range 0.03-0.05 we get results similar to TW, as expected 
due to the values of L. 

In the inner crust of a compact star, matter is not 
homogeneous, and is characterized by different isospin 
contents for each phase, i.e., the clusterized regions are 
more isospin symmetric than the surrounding nuclear 
gas, the so-called isospin distillation [28], l29(. The ex- 
tension of the distillation effect is model dependent and 
it has been shown that NL3 and other NLWM parameter- 
izations lead to larger distillation effects than the density 
dependent hadron models (30l - [32l |. In FigJ5]we show the 
ratio of the proton versus the neutron density fluctua- 
tions corresponding to the unstable mode. This ratio 
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defines the direction of the instability of the system. We 
show the results for two proton fractions Y p = 0.3, and 
0.05, for the sake of studying the effectiveness of the mod- 
els in restoring the symmetry in the liquid phase in two 
situations of interest for compact stars: /3-equilibrium 
matter with and without trapped neutrinos. The uj — p 
coupling decreases the Sp p /Sp n ratio, as already obtained 
in [22| for NL3cl>/9, although in this last work a dynam- 
ical calculation was performed. Comparing with the re- 
sult reported in [3(| we conclude that: a) QMC behaves 
differently from NLWM models such as NL3. For these 
models the ratio of the proton versus the neutron density 
fluctuations increases with the density, while for QMC af- 
ter a maximum obtained at p ~ 0.02 fm -3 , this ratio de- 
creases and more strongly if A„ is large; b) QMC presents 
a behavior similar to the one of relativistic models with 
density dependent couplings such as TW, however, the 
decrease of the distillation effect with density is not so 
strong [3(| HH , even for the largest A v coupling we have 
considered. 



B. Neutron stars 

Having discussed the behavior of the generalized QMC 
model at subsaturation densities, we now turn to the 
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FIG. 2. (Color online) Direction of instability for Yp=0.05 
(a), and 0.3 (b). 



main topic of the present work and discuss the stellar 
properties obtained in the present approach. The compo- 
sition of the stellar matter is determined by the require- 
ments of the charge neutrality and chemical equilibrium 
under the weak processes 

B 1 ->B 2 +l + v l ; B 2 + l^B 1 + v l (6) 

where B\ and B 2 are baryons, I is a lepton. The EOS 
which depends on the chemical potentials is now modi- 
fied according to [HJ , so that the lowest eight baryons are 
taken into account. As we restrict ourselves to zero tem- 
perature, no trapped neutrinos are considered, but the 
electrons and muons are considered so that charge neu- 
trality and /3-equilibrium can be enforced. The hyperon 
couplings are not relevant to the ground state properties 
of nuclear matter, but information about them can be 
available from the levels in hypernuclei [H, l3~4T - l3ci ] . Note 
that the s-quark is unaffected by the sigma and omega 
mesons i.e. g% = <j£ = . 

In QMC the couplings of the hyperons to the er-meson 
do not need to be hxed because the effective masses of 
the hyperons are determined self-consistently at the bag 
level. Only the x^b and x p b have to be fixed. The cou- 
pling strength of the p meson is given by the isospin of the 
baryon, and we obtain x u b from the hyperon potentials 
in nuclear matter, Ub = —(M B — Mb) + x u b9u^>o, for 
B = A, £ and 5 to be -28 MeV, 30 MeV and -18 MeV, 
respectively. We hnd that x^a = 0.743, i^s = 1.04 and 
%wE = 0.346. x p b = 1 is fixed for all the baryons. How- 
ever, while the binding of the A to symmetric nuclear 
matter is well settled experimentally [ijj], the binding 
values of the S~ and still have a lot of uncertain- 
ties [l6j]. We, therefore, test the effect of the coupling 
to the cascade and show results also for Vs = —10 and 
MeV. In fact, measurements from the production of 
5 in the 12 C(AT~, AT + )| 2 Be are compatible with a shal- 
low attractive potential Vs ~ — 14 MeV [l5| . We obtain 
i uH = 0.3989 for Vs = -10 MeV and x uB = 0.4643 for 
Vs = MeV. 

The resulting EOS are displayed in Figf3^,) for the 
QMC, QMCwp for different values of the coupling pa- 
rameter and QMC with protons and neutrons only. We 
also include the empirical EOS obtained by Steiner et 
al. from a heterogeneous set of seven neutron stars with 
well-determined distances [l9j]. We conclude that the 
agreement of the theoretical EOS with the empirical one 
when hyperons are included in the calculation is defined 
by the hyperon-meson interaction and the A v coupling, 
or, equivalently, by the symmetry energy. The QMC pn 
EOS agrees with the constraints. However, the inclu- 
sion of hyperons with the hyperon couplings obtained for 
the hyperon nuclear potentials taking Va = — 28 MeV, 
V s = 30 MeV and Vs = -18 MeV makes the EOS too 
soft. Increasing A v makes the EOS harder bringing the 
EOS closer to the constraints defined by the empirical 
EOS. This is easily understood with the help of FigJU 
Increasing A^ gives rise to a softer pn EOS at high den- 
sities and, therefore, hinders the onset of hyperons. So 
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FIG. 3. (Color online) a) Equation of state and b) TOV 
results for QMC and QMCwp. The EOS for V~ = -10 MeV 
and MeV where only obtained with A„ = 0.1. 



the larger the smaller the hyperon fraction in the star 
and the harder the EOS. 

The effect of a less attractive Vs potential is also clear: 
the EOS becomes harder because the onset of hyperons 
occurs at larger densities as shown in FigQ] We con- 
clude that any mechanism than hinders the formation of 
hyperons makes the EOS harder. 

The EOS enters as input to the Tolman-Volkoff- 
Oppenheimer [39j equations, which generate the macro- 
scopic stellar quantities. The obtained mass/radius curve 
for stars with a mass larger than 1M Q and the corre- 
sponding properties of maximum mass stars are then 
shown, respectively, in Figj3)D) and Table |TTJ 

First let us discuss the effect of the symmetry energy 
and the hyperon couplings on the mass/radius curve. A 
larger A„ gives rise to a softer EOS and, therefore, a 
smaller radius. It is seen that when going from A„ = 
to 0.1 the radius of a star with a mass M = 1 — 1.5 M 
decreases by ~ 0.3 Km. This effect was already discussed 
within NLWM for nucleonic stars [(| [23[ and for hyper- 
onic stars 24] . In this last paper it was shown that there 
exists a clear correlation between L and star radius. How- 



TABLE II. Stellar properties obtained with the QMC model 
and different values of the parameter A„ and the H-meson 
coupling, pn stands for nucleonic matter with no hyperons 
included. 
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FIG. 4. (Color online) Strangeness fraction as a function of 
density QMCc^p. For Vs = —10 and MeV we have taken 
A„ =- 0.1, For Vs = -18 MeV, we show results for A„ = 
0, 0.05, 0.1. 



ever, within the models discussed in [24J the maximum 
mass did not to depend on L, while in the framework of 
the QMC model there is a clear effect of almost O.IM© 
if A„ increases from to 0.1. This is mainly due to the 
smaller strangeness fraction inside the star. 

The reduction of the attractiveness of Vs has a similar 
effect on the maximum mass of the star, i.e., the mass 
increases ~ 0.2 Mq if Vh increases from -18 to MeV. 

We conclude that there is still quite a large uncertainty 
on the coupling of hyperons to nuclear matter and there- 
fore, there is still room for a very massive star such as 
the recently measured pulsar J1614-2230 with a mass 
M = 1.97 ± 0.04 [3], even including hyperons in the 
EOS. This, however, is a particularly massive star. Most 
of the known pulsars can be obtained by the present 
models. 



IV. SUMMARY AND DISCUSSION 

We have proposed a modified QMC model which in- 
cludes a nonlinear cop coupling in the same fashion as it 
has been proposed for the NLMW @. In QMC model 
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the nucleons are described as non-overlapping bags. The 
extra contribution allows the softening of the symmetry 
energy at large densities. In NLWM or QMC without this 
term the symmetry energy increases almost linearly with 
the baryonic density giving rise to very hard stellar mat- 
ter EOS. The inclusion of this term remedies this prob- 
lem and brings down the slope of the symmetry energy 
at saturation density to values closer to the experimental 
predictions (see Q for a compilation of all constraints 
on L). With the modified QMC, we may take advan- 
tage of the already known good properties of the QMC 
together with a symmetry energy that is not too hard. 
We have also shown that the behavior of the modified 
QMC model is closer to the properties of nuclear rela- 
tivistic models with density dependent couplings such as 
TW, a relatvisitic nuclear model with density dependent 
coulings [26j . Namely, the crust-core transition density is 
larger than the one predicted by the standard QMC and 
similar to TW and the distillation effect in non homoge- 
neous matter does not increase with density as in NL3, 
but decreases as in TW. 

We have also discussed the effect of the new EOS on 
the stellar properties. Hyperons were included in the 
EOS, and, for the hyperon couplings we took advan- 
tage of the fact that QMC predicts the hyperon effective 
masses without being necessary to fix the hyperon-u cou- 
plings. We have used information from hypernuclei to fix 
the hyperon-cj coupling and the hyperon- p coupling was 
considered equal to the one of the nucleon. Since there 



is a large uncertainty on hypernuclei with X and 3, we 
have also tested the effect of increasing the potential Vs 
so that it becomes less attractive. It was shown that both 
the symmetry energy and the hyperon couplings have a 
strong effect on the mass and radius of the star. A softer 
symmetry energy gives rise to smaller stars. Also the hy- 
peron fraction is affected: softer symmetry energy corre- 
sponds to a smaller hyperon fraction as already discussed 
in 24]. However, within QMC the density dependence of 
the symmetry energy has also an effect on the maximum 
star mass and this effect was not obtained in 24] . 

It was also shown that the hyperon nuclear interac- 
tion defines the amount of strangeness in the star, and, 
therefore, has a strong influence on the maximum mass 
allowed. Even including hyperons in the QMC EOS we 
could explain the mass of the pulsar J1614-2230 if the 
cascade nuclear potential is set to be very little attrac- 
tive. More data on hypernuclei is needed to constrain the 
hyperon-meson couplings. 
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